Empirical transverse charge densities in the deuteron 
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Using the recent empirical information on the deuteron electromagnetic form factors we map out 
the transverse charge density in the deuteron as viewed from a light front moving towards the 
deuteron. The charge densities for a transversely polarized deuteron are characterized by monopole, 
dipole and quadrupole patterns. 



OO 

o 
o 

(N 



00 
(N 



Oh 
I 

Oh 

CD 



> 

en 
vn 
^, 

o 

00 

o 



X 



PACS numbers: 13.40.Gp, 21.10.Ft, 21.45.Bc 

Electromagnetic form factors (FFs) of the deuteron 
have received a lot of attention in recent years; for recent 
reviews see e.g. Refs. [lll3]- Based on the large amount of 
precise deuteron form factor data, it is therefore of inter- 
est to exhibit the spatial information on the quark charge 
distributions in the deuteron. In this letter, we develop 
the general formalism to extract charge densities for a 
spin-1 particle and apply it to the case of the deuteron. 

In the following we consider the electromagnetic 
deuteron elastic FFs when viewed from a light front 
moving towards the deuteron. Equivalently, this cor- 
responds with a frame where the deuterons have a 
large momentum-component along the z-axis chosen 
along the direction of P = (p + p')/2, where p {p') 
are the intial (final) deuteron four-momenta. We indi- 
cate the deuteron light-front + component by P+ (defin- 
ing fl^ = a^ ± a^). We furthermore choose a light-front 
frame where the virtual photon four-momentum q has 
q^ — 0, and has a transverse component (lying in the 
xy-plane) indicated by the transverse vector cf^, satsi- 
fying q^ = ~'ti_ — ~Q^- Ir* such a light-front frame, 
the virtual photon only couples to forward moving par- 
tons and the + component of the electromagnetic cur- 
rent /+ has the interpretation of the quark charge den- 
sity operator. It is given by J"'" (0) = +2/3m(0)7+m(0) — 
l/3d(0)7+(0)d(0), considering only u and d quarks. 
Each term in the expression is a positive operator since 

In the following, we will use empirical information on 
deuteron elastic FFs to study the deuteron quark charge 
densities in the transverse plane. It is customary to 
denote the three deuteron elastic e.m. FFs by Gq, Gm, 
and Gq, corresponding to the Coulomb monopole (Gc), 
magnetic dipole (Gjvi), and Coulomb quadrupole (Gq) 
FFs respectively. Similar relations between nucleon den- 
sities and FFs can be found in [3, 4, 5, 6]. 

We start by expressing the matrix elements of the 
/+(0) operator between deuteron states as. 



where A = ±1, (A' = ±1,0) denotes the initial (fi- 
nal) deuteron light-front helicity, and where q±^ = 
Q{cos^q Sx + siiKpq Cy). Furthermore in Eq. ^, the he- 
licity form factors Gt, , are real (due to time reversal in- 
variance) and depend only on Q^. 

We can then define transverse charge densities for 
deuteron in helicity states of A = ±1 or A = as. 
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where b = b (cos (pi, e^- + sin </)(, iy) denotes the position 
in the xy plane from the transverse cm. of the deuteron. 
The two independent helicity conserving FFs Gjj^j and 
Gqq can be expressed in terms of Gc,m,Q ^s. 
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with t] = Q^/(4M^), and M^ is the deuteron mass. 

The definitions and normalizations of Gq, Gm, and 
Gq are the customary ones obtained from the matrix el- 
ements of the electromagnetic current 17|], 
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with £i 2 the polarization vectors of the initial and final 
deuteron. The charge and quadrupole FFs follow from. 



Gc = Gi + -Tl Gq, 



Gi + (1 + rj) G3 - Gm, 
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with normalizations Gc(0) = 1, Gm(0) = }id [mag- 
netic moment in units e/(2Mrf)], and Gq(0) = Q^ 
(quadrupole moment in units e/M^). 



For numerical evaluation, we use the parameteriza- 
tion of the deuteron form factor data given as fit II by 
Abbott et al. |3]. This parameterization is based on forms 
suggested in H, which read 
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where the dipole form factor G{Q^) = (1 + Q^/S^)~-^ 
has a (non-standard) mass parameter S = 898.52 MeV. 
The reduced amplitudes are 
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and the values of the parameters obtained by the au- 
thors of 1 8] are given in Table ID 
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TABLE I: Parameters for deuteron form factor fit II of Ref. [g]. 

Figure[T]shows the transverse charge densities for def- 
inite helicity (longitudinally polarized) deuterons. Re- 
call that the transverse charge density is the charge den- 
sity projected onto a plane perpendicular to the line of 
sight, which also defines the longitudinal direction. The 
charge density for the A = 1 state, in the upper panel, 
is smooth and peaks in the center. The A = state, in 
the middle panel, features a dip in the center, which can 
also be seen in the lower panel, where the density along 
the y-axis is plotted for both states. 

This dip — found directly from the data — reflects the 
low charge density or "hole" in the center of the 
deuteron wave function. For the helicity-0 state, calcu- 
lated S- and D- state wave functions produce toroidal 
equidensity surfaces in the central regions, with the axis 
along the z direction jlCl||. Hence looking along the z-axis 
for this state, one can see the hole in the center, albeit the 
hole be partly filled in because one must look though 
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FIG. 1: Quark transverse charge densities in the deuteron. Up- 
per panel : pf. Middle panel : pjj. The lightest (darkest) regions 
correspond with largest (smallest) densities. Lower panel : the 
density along the y-axis for p'^ (dashed curve) and pj] (solid 
curve). For the deuteron FFs, we use the empirical parameter- 
ization II of Ref. 1 8], based on Ref. ^. 



the outer layers of the deuteron which are mostly S-state 
and spherically symmetric. We emphasize that the dip 
seen in our figure is based only on observation (codified 
in the form factor fits) and a light-front interpretation 
of the data. The result is consistent with but does not 
use calculations based on models of nucleon-nucleon 
interactions. Further, the light front viewpoint gives a 



two-dimensional charge density that is relativistically 
correct, unlike charge densities obtained from Fourier 
transforming form factors in an equal time formalism. 

We next consider the charge densities for a trans- 
versely polarized deuteron, denoting the transverse po- 
larization direction by Sj^ = coscpg Cx + sin ^5 Cy. The 
transverse charge densities can be defined as. 
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where Sx is the deuteron spin projection along the direc- 
tion of Sx • The transverse spin state can be expressed in 
terms of the light front helicity spinor states as, 

+ e2%|A= -1)), 



|sx=0) = -^(|A = +l>-e2'^MA = -l>).(9) 

By working out the Fourier transform in Eq. lO for the 
two cases where sx = +1 and sx = 0, one obtains, 

/ti(^) = ^^Q(/o(^'Q)^(G++G+^ 



+ s\n{^i-(Ps)h{hQ)V2G+^ (10) 
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+ cos2{(p^-cPs)WQ)Gt^^^). (11) 

One notices from Eqs. lllOlllll that the transverse charge 
density pj-^ contains monopole, dipole and quadrupole 
field patterns, whereas Pj^ only contains monopole and 
quadrupole field patterns. The deuteron helicity FF with 
one unit of helicity flip, which governs the dipole field 
pattern in jO^ ^ , can be expressed in terms of Gq^m.q as, 

^01 = - YT^ {gc - ^(1 - 7) Gm + I Gq}, (12) 

whereas the deuteron helicity FF with two units of he- 
licity flip, corresponding with the quadrupole field pat- 
terns, can be expressed as. 
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Note that the four deuteron helicity FFs introduced here 
are not independent. The angular condition jllll relating 
them reads, 

(2^ + 1)G+ + 2^2^ Go+ + G+i +1 - G+, = 0. (14) 



It is instructive to evaluate the electric dipole moment 
(EDM) corresponding with the transverse charge densi- 
ties jOjg , which is defined as. 
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Eq. ||TTJ gives dg 
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Expressing the spin-1 magnetic moment in terms of the 
^-factor, i.e. Gm(0) = g, one sees that the induced EDM 
dj is proportional to g — 2. The same result was found 
for the case of a spin-1 /2 particle [5]. One thus ob- 
serves that for a particle without internal structure (cor- 
responding with g = 2 [12]), there is no induced EDM. 

One can also evaluate the electric quadrupole mo- 
ment corresponding with the transverse charge densi- 
ties pj^ . Choosing Sx = Sx, the electric quadrupole 
moment can be defined as. 
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From Eqs. l llOllll l one obtains, 
1 



Qi 



tQo 



= l{[GM(0)-2] + [GQ(0) + l]}f^j.(18) 

We may note that for a spin-1 particle without inter- 
nal structure, such as the W and Z gauge bosons in the 
standard electroweak theory, it is required that at tree 
level Gm(0) = 2 and Gq(0) = —1, in order to satisfy 
the Gerasimov-Drell-Hearn sum rule to lowest order in 
perturbation theory [13, 14]. For the gauge bosons, any 
deviations from these values would indicate new (be- 
yond standard model) physics. Current data from pp 
collisions (up to y/s = 1.96 TeV) are compatible with 
zero anomalous WWj and ZZ7 couplings [15]. 

It is thus interesting to observe from Eq. JTSl l that Q^ 
is only sensitive to the anomalous parts of the sprn- 
1 magnetic dipole and electric quadrupole moments, 
and vanishes for a particle without internal structure. 
The deuteron has a magnetic dipole moment Ga4(0) = 
1.71 [16], close to the natural value for a spin-1 particle. 
However, in contrast to the W and Z gauge bosons, the 
deuteron has a large anomalous quadrupole moment. 
Its measured value is Gq(0) = 25.84 (3) [17], highlight- 
ing the prominent role of the pion exchange potential. 
Eq. | [T8)l thus shows that for a deuteron polarized along 
the X-axis, its transverse charge densities show large 
quadrupole moments (Qj > and Qg < 0). 

Pictorial results for the transverse charge density with 
transverse deuteron polarization are shown in Fig. [J] 
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FIG. 2: Quark transverse charge densities for a deuteron which 
is polarized along the positive x-axis. Upper panel : pj-y. Mid- 
dle panel : Pj„. The light (dark) regions correspond with 
largest (smallest) densities. The lower panel compares the den- 
sity along the y-axis for pj^ (dashed curve) and p^g (solid 
curve). For the deuteron e.m. FFs, we use the empirical pa- 
rameterization 11 of Ref. LSI, based on Ref. l^. 



The upper panel shows the charge distribution for the 
state with projection-1 in the x-direction. One sees the 
large effects of the quadrupole term together with a 
small overall shift of the charge distribution alon g the y- 
axis. The latter is consonant with the observation il8lll9ll 
that an object with a magnetic dipole moment when sta- 



tionary exhibits an electric dipole moment when mov- 
ing, or when seen by a moving observer, proportional to 
the vector product v x (magnetic moment). The middle 
panel shows charge density for a state with projection- 
in the x-direction. The quadrupole term stretches the 
charge along the y-axis but does not cause any shift of 
the charge center. The two cases are compared in the 
lower panel, which plots the density along the y-axis. 

In summary, we used recent empirical information on 
the deuteron e.m. FFs to map out the transverse charge 
densities in longitudinally and transversely polarized 
deuterons. Notably, one sees a dip in the center of the 
charge distribution for helicity-zero deuterons. This is in 
harmony with nuclear force model calculations, which 
for the zero helicity case predict toroidal equidensity 
surfaces for higher densities, with axis along the quan- 
tization direction. Transversely polarized deuterons 
show dipole and quadrupole structure in the charge dis- 
tributions. Their electric dipole and quadrupole mo- 
ments only depend on the spin-1 particle's anomalous 
magnetic dipole moment and its anomalous electric 
quadrupole moment, arising from its internal structure. 
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